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We study the competition between pinning of a charge density wave (CDW) by random dis- 
tributed impurities and a periodic potential of the underlying crystal lattice. In d = 3 dimensions, 
we find for commensurate phases of order p > pc ~ O/tt a disorder driven continuous roughening 
transition from a 'flat' phase with translational long range order to a 'rough' glassy phase with quasi 
long range order. Critical exponents are calculated in a double expansion in ^ — /pi — 1 and 
e = 4 — d and fulfill the scaling relations of random field models. Implications for flux line lattices 
in high temperature superconductors are briefly discussed. 
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Charge density waves (CDWs) are broken symmetry 
states of mostly strongly anisotropic quasi-one dimen- 
sional metals Their ground state is a coherent 
superposition of electron-hole pairs and shows a peri- 
odic spatial modulation with wave vector qo — 2\<ip. 
The Fermi wavevector k^? enters here due to the diver- 
gent electronic response at 2kj^ in one dimension. This 
in combination with a sufficiently large elect ron-phonon 
coupling may lead to a Peierls instability, resulting in a 
new state with a periodically varying charge modulation 
of density 



p(x) « po + Pi cos(qox -f 0(x)). 



(1) 



For a partially filled electron band, the period 7r/fci? is 
in general incommensurate with the underlying lattice. 
Considering the excitations of CDWs well below the in- 
stabihty temperature, one may neglect variations of the 
amplitude pi since they are accompanied by relatively 
high energy expenses. In contrast, long wavelength defor- 
mations of the phase (/'(x) are accompanied by arbitrarily 
small elastic energy costs. In perfectly clean crystals the 
phase excitations form a gapless excitation branch, giving 
rise to an infinite conductivity at zero frequency. 

A very similar picture is found for other incommen- 
surate structures like spin density waves in anisotropic 
metals [Q, mass density waves in superionic conductors 
1^, polarization density waves in incommensurate ferro- 
electrics or fiux-line arrays in type-II superconductors 
1^. In the absence of a coupling to the underlying lattice, 
the wavevector qo of the modulation is related to external 
parameters like temperature for ferroelectrics or the ex- 
ternal magnetic field for superconductors and in general 
incommensurate with one of the reciprocal lattice vectors 
g of the crystal. Although our considerations apply to 
this large group of systems as well, we will continue to 
use the terminology of CDWs. Its extension to the other 
cases is straightforward and will be discussed for fiux- 
line lattices in high-temperature superconductors briefiy 
at the end. 

Lattice eflFects, which have been neglected so far, may 



become relevant if qo is close to g/p where p is an integer 
corresponding to 1/p filling of the conduction band. For 
a sufficiently strong lattice potential g^v (in units of the 
phase stiffness constant 7 ), the density wave locks in at 
the commensurate wavevector g/p. In this case phase 
ffuctuations are suppresed, i.e. the CDW is in a 'fiat' 
phase. The phason branch acquires a gap of order v and 
an applied external force density /cx has to overcome a 
threshold value fc,v ~ ig^pv to depin the structure. If g 
and qo are parallel, the width of the commensurate phase 
is given by \g/p - go| = i5 < « gv'^/'^. 

The actual modulation of ^(x) in these phases is in 
general very complicated. It is therefore convenient to 
go over to a new phase field which describes the distor- 
tions from the lock-in wavevector, such that (j) — Q in 
the commensurate (C) phase. The incommensurate (I) 
phase close to the lock-in is then described by a lattice of 
domain walls of intrinsic width ^0 ~ '^1 {pgv^^"^) and di- 
verging separation I if the C phase is approached. There 
is no lattice pinning in the I phase. 

Thermal fiuctuations strongly diminish C phases in 
two dimensions In three dimensions thermal fiuctu- 
ations have only minor effects close to the transition to 
the disordered phase. On the other hand, experimentally, 
commensurability effects have rarely been seen, which 
was mainly explained by the smallness of v. 

Disorder is another ingredient which is known to 
change this simple picture. If the potential of an indi- 
vidual impurity is weak, it is sufficient to consider its 
effect on the phase order only. Neglecting commensura- 
bility effects, it was shown by Fukuyama and Lee ||^ that 
long range order of phases (/)(x) is indeed destroyed by an 
arbitrarily weak random potential in dimensions d < 4 
on length scales L larger than the Fukuyama-Lee-length 
La = A-i/(4-d)^ ^ ^ -i?"(0) denotes the strength 
of the disorder fluctuations, where R{(j)) is defined in Eq. 
(^) . As a result of the adaption of the CDW to the im- 
purities, the modulated structure is now also pinned in 
the I phase with a pinning force density /c,a ~ 7^a^- 

In general, however, lattice and impurity pinning will 
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occur simultaneously, but act in opposite directions. 
Whereas the lattice prefers registry of the CDW, result- 
ing in a commensurate flat phase, the disorder tends to 
destroy its long range order and makes the CDW rough. 
It is therefore the aim of the present paper to study the 
combined effect of both pinning mechanisms. We will re- 
strict ourselves in the following to the realistic situation 
where gL\,gS^Q 3> 1, i.e. to the weak pinning case. Since 
pinning forces from the lattice and disorder act in gen- 
eral in opposite directions with maximal values fc.v and 
/c,A, respectively, one naively expects that the C phase 
is stable if /,,A < /c,., i.e. for A < (g»(4-'^)/2 |]. 

A more detailed calculation to be presented below how- 
ever shows that for 5 — Q andp < Pc{d) (pc(3) « 6/7r) the 
pinning due to the lattice is dominating even for arbitrar- 
ily small potential values v, and the C phase persists up 
to an intermediate disorder strength. Inside the C phase, 
the translational order parameter ^'(x) = pi exp(i0(x)) 
takes a non-zero expectation value (^'(x)), where ( ) and 
the bar denote thermal and configurational average, re- 
spectively. For p > Pc{d), the system indeed undergoes 
a disorder driven continuous roughening transition if A 
exceeds a threshold value Ac ^ v'^^^ with e = 4 — d. For 
A < Ac, |(^'(x))| (Ac - A)'^, and the correlation func- 
tion is:(x) = ((5*(x)^**(0)), (5^'(x) 'J'(x) - (*(x)), 
shows an exponential decay. In contrast, for A > Ac 
the phase is rough and unlocked despite ^ = with 
(\I'(x)) = but a power law decay of K{x.) ^ |x|^~'^~^. 
From a functional RG calculation we find in a double ex- 
pansion Q in e <C 1 and fi — p^ /p^ — I <ti 1 the critical 
exponents 



= 2-e, 



(2) 



where v denotes the correlation length exponent in the 
C phase. The exponent 9 is the negative temperature 
eigenvalue which appears in the random field hyperscal- 
ing relation v{d — 6) = 2 — a [|o|. The critical exponents 
fulfill all exponent relations and inequalities known from 
random field models with discrete symmetry of the order 
parameter [nOl. In particular, we obtain the exponent 
f] = e(i + 7r279) at Ac. 

Our starting point for the description of the CDW close 
to a C phase with a p-fold degenerate groundstate is the 
following Ginzburg-Landau-Hamiltonian H 



n=^jd''^\^]^{V4>-8f 



g\coB{p4,) + V{<j),^)\ (3) 



with 6 = 6z. Here the first term represents the elas- 
tic energy of the long-wavelength phase deformations 
controlled by the stiffness constant 7 = n{eF)h'^Vp/2 
with Fermi velocity v-p and density of states n{e-p) at 
the Fermi level, g^v ^ p\ ^ measures the strength 
of the lattice potential relative to 7. The last term 
F((/!),x) = piVbY^cos(0(x) — a(x)) describes the effect 
of impurities, where 7V0 and Ui are the strength and 
the density of impurities, respectively. The uniformly 
distributed and uncorrelated random phase a(x) arises 



from the random impurity positions since we have as- 
sumed the weak pinning case with n^L^ ^ 1 The 
disorder correlation function R{(j)) is given by 



V{(j), x)T^(0', x') = - (/)')(5(x - x') 



(4) 



with a bare function RQ{(f>) — nip\VQ cos((/))/2. 

We first estimate the influence of the pinning poten- 
tials by perturbation theory. For a vanishing lattice po- 
tential u = 0, to begin with, it has bee n shown [pT| that 
the phase correlation function C(x) — {{(f>{x) — 0(0))^) 
behaves logarithmically on length scales |x| > La, i.e. 
C(x) = 2A(4- (i)ln|x/LA| with A = w 1.1 [||. 

Neglecting the non-Gaussian character of which is 

justified for e ^ 1, lowest order perturbation theory in v 
yields an effective Hamiltonian with mass |13 



g^p^ vcos} 



-p^A(4-d)/2 



(5) 



on length scales L > La. Comparing this power of L 
with the L^^ behavior of the elastic term, we conclude 
that the perturbation is always relevant if p < Pc = 
(4/A(4- For lar ger p a small periodic potential 

can be neglected. 

Next we employ an RG calculation, which starts with 
the disorder averaged replica-Hamiltonian 



_L 

2T 



af3 



(6) 



A^ = {Tr/3)'^eA''K~^, which describes 



In the absence of crystal pinning, v = 0, the misfit S can 
be shifted away by the transformation ^(x) ^(x) -I- Sz 
due to the statistical symmetry of the disorder under 
spatial translations. The corresponding Hamiltonian has 
been studied repeatedly in the past (l^,|l^ using the func- 
tional RG-technique in c? = 4 — e dimensions. The result- 
ing RG flow leads to a zero-temperature fixed point U* 
where -i?^(0) 

the glassy unlocked phase (U) of the CDW. 

To extend the calculation to w 7^ 0, we follow here the 
functional RG treatment developed by Balents and Kar- 
dar for a closely related case |l^ . Although for p > pc a 
weak periodic potential is irrelevant, one expects that it 
becomes a relevant perturbation once v exceeds a critical 
value Vc. Since Vc has to vanish if p —^ p'^ , it is tempting 
to describe the new fixed point F* corresponding to this 
transition perturbatively in e = 4 — c? and /i = p'^ /Pc ^ 1- 
A finite value of v close to F* will lead in particular 
to a renormalization of 7/T. While keeping 7 fixed, we 
absorb this in an additional renormalization of T, R{4') 
and V. The fixed point value /S.*p of A at F* differs from 
A^ by corrections of order e/i. The renormalization of 
T will hence not change the sign of —6* = —2 + 0{e,fi), 
describing the irrelevance of T also at F* . 
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(a) A 




FIG. 1. (a) Schematic RG flow for 5 = and p > pc in 
the u-A-plane. A critical hne Ac{v) with fixed point F* sep- 
arates the flat ((7) and rough, unlocked ( U) phase. For p < pc 
the fixed points F* and U* merge. Since for large disorder 
strength A topological defects will proliferate, we expect a 
disordered phase D above the line Ad{v) but the RG flow in 
this range is not yet clear. The transition from the C to the D 
phase is probably in the universality class of the random field 
p-state clock model, (b) Possible phase diagrams for finite S 
at fixed V for p > pc and (c) for p < pc- 

Under the RG transformation higher order harmon- 
ics VmC0S.{pm4>), m > 1, will be generated. However, 
the Vm with m > 1 are irrelevant at F* since the eigen- 
values A„i = 2 — ^p'^m'^Ap = 2(1 — m 
these pertubations are negative for fi 
We can restrict ourselves therefore to the lowest order 
harmonic vcos{p(j)). Note that also terms of the form 
{R"{(t)a — cos{mp(j)a), generated by the RG, are 

strongly irrelevant with scaling dimension A„i — d — 2. 
Collecting the contributions to the RG-flow up to the 
first non- linear order in v and R{(t>), we obtain the fol- 
lowing RG flow equations valid close to the T = fixed 
points and to lowest order in e and /x with A = — i?"(0). 



2a + of 
if m > 1. 



dT 

m 

dv 
'dt 
dR{(f>) 
dt 



^ = i2-d-—Av 



2-P-A 
2 

„6 



{v,A)T 



^A.3 



(7) 
(8) 



e-^z;2A|i?(0) + ii?"(0f -f Ai?"(0) (9) 



substitutions A '^g'^v v, 
i^-i = 2''-i7r'^/2r(d/2). Next 



where we made the 
KdA'^-'^Ri^) ^ i?(</)), 
we analyze the fixed points of the RG flow. Rewrit- 
ing R{<p) in a Fourier series — '^'^^ 
the fixed point function RIj{4>) ||12[ for w = (as well 
for any other fixed value of v) can be parameterized by 
^m^Rl^ = A*(t;). Rewriting i?™ = -^m-'^A + r„, the 
are irrelevant with eigenvalues — e(l — ii)cm with pos- 
itive factors Cm ^ tti^- Since we are interested in the 
universal properties of the transition we will therefore 
adopt here the particularly simple functional form of the 



A^. 



(10) 



fixed points which holds approximatively on length scales 
L > La- Eq- @ reduces then to a single flow equation, 

-;- = ^A^ - '-- 

dt 7r2 2 

For finite v, the resulting phase diagram is character- 
ized by a critical line separating the flat and rough 
phases, along which the RG flow is driven into a new 
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critical T — fixed point with A|^. — ^e(l — ^) and 
Vp — ^J7\ij2jp^ (see Fig. |^). Linearizing around this 
fixed point, we obtain the critical exponents of Eq. (|^). 
This transition can be approached e.g. by changing the 
temperature T since v^l'^ j A ~ p}"*^^^^ where the CDW 
amplitude p\ changes with T and vanishes at the Peierls 
temperature. 

So far we have considered the case (5 = 0. For non-zero 
(5 the linear gradient term in (^) can be integrated out 
leading to —8^ J d'^~'^yA(j){y) where A0(y) = (/)(y, L) — 
(j){y, 0). In the C phase A(j){y) = and hence the 6 — 
calculation applies. To estimate dc, we compare in (^) 
the term 5d<j)/dz with the cosine term on the scale which 
corresponds to the renormalized soltion width (or corre- 
lation length) ^ leading to 5c « IM = gv^/'^{l- A/ A^f . 
In the I phase we rewrite A(f) — 2^-^ and treat I as a pa- 
rameter to be determined later by minimizing the free 
energy. With the transformation (/)(x) — 0(x) — we 
go back to periodic boundary conditions [ p^ changing the 
cosine term to g'^v cos {pcf) — ^j^). If ^ ^ ^ the RG up to ^ 
is insensitive to the oszillations in z and the soliton lattice 
can be described in terms of an effective free energy which 
includes also the steric repulsion between the walls due 
to random field like disorder. The resulting domain wall 
density vanishes as ^ {5 — Sc)^ by approaching the C 
phase where ^ = C/2(l-C) = (5-d)/2(d-2) Jll follows 
from the roughness exponent Q = {5 — d)/3 |l^ . More- 
over, the random potential destroys the translational long 
range order of the soliton lattice, which shows an alge- 
braic decay of correlations beyond ^ I and ^|| ~ 



if(x)=p?cos^| (^) +(- 




-C7r^/18p^ 



(11) 



In the opposite case ^ I, v averages to zero on scales 
larger than I for geometric reasons. The correlations 
show a crossover to a behavior which is given by ( |ll| ) 
with ^_L, — > La and exponent — ett^/IS. 

As has been discussed earlier 20 , a small applied 
external force of density /ox = eE from an electric field 
applied along the chain direction leads for /ex /c, to a 
creep motion of the CDW with a creep velocity 



v{fex) ~ exp 



El 

' T 



(12) 



2 and 



In the U phase, k = (d— 2)/2, /^ = /c,a and E, 
whereas in the C phase k = d — 1, /c = fc,v ~ 7^ 
Ec ~ 7^^^^. Because of the pronounced difference in the 
exponents k, measurements of the creep velocity should 
give clear indications about which phase is present. 
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Recently measured I-V curves of the conductor o-TaSs 
at temperatures below IK can be fitted by (^2|) with 
K — 1.5 - 2 The experimentally observed tendency 
to larger k for purer crystals confirms the above inter- 
pretation. In several materials, such as K0.3M0O3, the 
periods are near p — A commensurability at low tem- 
peratures. For this material one obtains ^0 ~ 10~^cm 
l^. The typical parameters jVo « 10~^eV, pi = 10~^ 
and lO^cm/sec yield, after proper rescaling of 

anisotropy, the estimate La ~ 10~^cm for an impurity 
density of lOOppm. Thus it should be possible to see 
commensurability effects if S becomes small enough. 

So far, we have excluded topological defects. These 
can be considered if we treat 0(x) as a multivalued field 
which may jump by ±2tt at surfaces along which electrons 
are inserted to / omitted from the chains at the intersec- 
tion points of the chains and surfaces. These surfaces 
are bounded by vortex lines. For S ^ = v it has been 
argued recently that for weak enough disorder strength 
A < Ad, the system is stable with respect to the forma- 
tion of vortices However, vortex lines will prolif- 
erate for A > Ajj. At present it is not clear whether the 
corresponding transition is continuous or first order. For 
S — but V > we expect that this transition extends 
to a line Ad{v) until v reaches a critical value vd with 
Ad{vd) = Ac{vd) (see Fig. |l|). For larger v the tran- 
sition is probably in the universality class of the p-state 
clock model in a random field, which has an upper critical 
dimension dc = 6. In the Ising casep = 2 the transition is 
of second order . A non-zero value of S will in general 
increase the size of the I phase, as schematically scetched 
in Fig. |l|b,c. 

Our results can also be applied to the pinning of flux 
line lattices (FLL) in type-II superconductors. In lay- 
ered superconductors, the Cu02 planes provide a strong 
pinning potential favoring, for flux lines oriented parallel 
to the layers, a smectic phase with translational order 
present only along the layering axis 1^ . The influence of 
disorder on this phase is described by the CDW model 
studied in this paper, if the CDW phase 0(x) is identi- 
fied with the deviations of the smectic layers from their 
locked-in state. Also a FLL oriented perpendicular to the 
layers in general feels a very weak periodic potential of 
the underlying crystal, but now the flux line displace- 
ments are described by a vector field. Since also in this 
case weak disorder leads only to logarithmically grow- 
ing displacement of the flux lines [p^ , a disorder driven 
roughening transition of the CDW type studied above 
can be expected for the FLL. Furthermore, numerical 
simulations of FLL's are usually modeled on a discrete 
lattice which acts for the FLL like a periodic pinning po- 
tential. Thus the used grid size has to be small enough 
(sufficiently large p) to avoid commensurability effects. 

In conclusion, we have shown that a 3D CDW which 
is subject both to impurity and lattice pinning under- 
goes a roughening transition from a flat commensurate 
to a rough and in general incommensurate phase. On 
the C side, this transition is described by three indepen- 
dent critical exponents which obey the scaling relations 



for random field systems with a discrete symmetry of the 
order parameter. Compared to the thermal roughening 
transition in = 2, where the T axis represents a line 
of fixed points, in the disorder driven case the existence 
of one universal fixed point leads to a different RG 
flow. Notice that in the disorder dominated system stud- 
ied here, the important logarithmic roughness sets in only 
beyond the Fukuyama-Lee length La, whereas thermal 
roughness in d = 2 appears on all length scales. The 
changing creep behavior should be an interesting conse- 
quence for the experimental observation of the transition. 
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